Thermal escape of fractional vortices in long Josephson junctions 
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We consider a fractional Josephson vortex in a long 0-k Josephson junction. A uniformly applied 
bias current exerts a Lorentz force on the vortex. If the bias current exceeds the critical current, 
an integer fluxon is torn off the K-vortex and the junction switches to the voltage state. In the 
presence of thermal fluctuations the escape process takes place with finite probability already at 
subcritical values of the bias current. We experimentally investigate the thermally induced escape 
of a fractional vortex by high resolution measurements of the critical current as a function of the 
topological charge k of the vortex and compare the results to numerical simulations for finite junction 
lengths and to theoretical predictions for infinite junction lengths. To study the effect caused by 
the junction geometry we compare the vortex escape in annular and linear junctions. 

PACS numbers; 05.45.-a, 74.50.+r 85.25.Cp 
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I. INTRODUCTION 

Vortices in long Josephson junctions (LJJs) usually 
carry a single magnetic flux quantum $0 and therefore 
are called fluxons. The study of fluxons has been attract- 
ing a lot of attention during the last 40 years because of 
their interesting nonlinear nature [ll [3, y] as well as be- 
cause of possible applications [1, H, H, 131 ■ 

Recently it was shown that one can create and study 
experimentally vortices that carry only a fraction of the 
magnetic flux quantum. Initially, vortices carrying only 
half of the magnetic flux quantum, so-called semifluxons 
d, [^, were observed and studied. They exist in 0- 
TT LJJs, which can be fabricated using superconductors 
with an anisotropic order parameter that changes sign 
depending on the direction in k space (e.g., d-wave order 
parameter symmetry [ll|, [l^, [H, 13 ) or with an oscillat- 
ing order parameter (e.g., with a ferromagnetic barrier 

In O-TT LJJs, some parts behave as junctions and other 
parts as tt junctions. In this kind of junctions, the ground 
state phase fj,{x) will have the value deep inside the 0- 
region and the value tt deep inside the 7r-region. In the 
Aj-vicinity of the O-tt boundary, fi{x) will continuously 
change from to tt, where Aj is the Josephson pene- 
tration depth. Because of the phase bending a magnetic 
fleld oc djji/dx localized in the Aj-vicinity of the boundary 
appears and a supercurrent (x sin[/i(a;)] circulates around 
it (1,1 The total magnetic flux localized at the O-tt 
boundary is equal to -|-<I>o/2. Since the Josephson phase 
is 27r-periodic, its value can be — tt inside the 7r-region 
instead of tt. In this case, the supercurrent flows counter- 
clockwise and an antisemifluxon with the magnetic flux 
—$0/2 is formed. 

It turns out that instead of a 7r-discontinuity of the 
Josephson phase at the O-tt boundary one can artificially 
create any arbitrary K-discontinuity at any point of the 
LJJ and the value of k can be tuned electronically p^ . 
At the resulting 0-k boundary two types of vortices may 



exist [T3|: a direct vortex with topological charge —k 
{—K vortex) and a complementary vortex with topologi- 
cal charge 27r — k {2tt — k vortex) . Without loosing gen- 
erality we assume that < k < 2tt. Such vortices are 
a generalization of semifluxons and antisemifluxons dis- 
cussed above. Only if their topological charge is less or 
equal than 2n by absolute value, they are stable [isj . 

By applying a spatially uniform bias current to a LJJ 
with a K vortex, one exerts a Lorentz force which pushes 
the vortex along the junction. Depending on the mutual 
polarity of the vortex and the bias current, the direction 
of the force can be adjusted. Since the vortex exists to 
compensate the phase discontinuity, it is pinned at it, 
i.e., it may only bend under the action of the Lorentz 
force, but it cannot move away from the discontinuity. 
However, increasing the bias current, the force can be 
made strong enough to tear off a whole integer fluxon out 
of a K vortex. The fluxon moves away along the junction, 
leaving a complementary [k — 27r) vortex pinned at the 
discontinuity. Further dynamics leads to the switching 
of the 0-K LJJ into the voltage state. In annular LJJs, 
this process takes place when the normalized bias current 
7 ~ I/Ico reaches the critical current of p^.[20| 



7c (k) 



sin(K/2) 



k/2 



(1) 



where Ico ~ jcwL is the "intrinsic" critical current, which 
corresponds to the measurable critical current ii k = 0; w 
is the LJJ width and L is its length. Due to thermal fluc- 
tuations, the escape process described above will already 
take place with finite probability at 7 < 7c. In this pa- 
per, we study the thermal escape of an arbitrary k vortex 
experimentally and compare it to numerical simulations 
for finite junction lengths and to theoretical predictions 
for infinite lengths. 

The paper is organized as follows. In Sec.|TTl we briefly 
review the thermal escape of the phase in a point-like 
Josephson junction. Then, in Sec. Illlj we discuss theo- 
retically the phase escape in a long JJ with a vortex and 
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compare it with our experimental findings in Sec. IIVI 
The influence on junction geometry is further investi- 
gated in Sec. |Vl 



II. POINT-LIKE JOSEPHSON JUNCTION 



In the Stewart-McCumber model [2l|,|22|, see Fig.IHa), 
the dynamics of a current-biased point-like JJ is de- 
scribed by an equation of motion for the Josephson 
phase 



-ill" 



1 

R\2Tr 



fi + Ic sm{fj.) 



0, (2) 



where C is the effective shunt capacitance, R is the ef- 
fective shunt resistance of the JJ embedded in the bias 
circuitry and dots denote the derivative with respect to 
time. The equation of motion ([2]) is equivalent to the 
damped motion of a particle of mass m = C($o/27r) 
along the generalized coordinate /i in a tilted washboard 
potential 



U{n) = -£',/(7/i + cos/i), 



(3) 



where Ej ~ $o/c/27r is the Josephson coupling energy, 
see Fig. [TJb). The tilting angle of the potential is pro- 
portional to the bias current 7 = I / Ic applied. 



a) 



1 
J 




FIG. 1: (a) Resistively and capacitively shunted junction 
(RCSJ) model, (b) Particle in a tilted washboard potential. 
The thermal escape from the metastable state is indicated by 
an arrow. 

The process of the particle escape from a metastable 
minimum in this Kramers-like system j23l [23 | corre- 
sponds to a transition of the JJ from a superconduct- 
ing zero-voltage state to a finite voltage state. In the 
absence of thermal and quantum fluctuations the junc- 
tion switches from the zero-voltage state, which corre- 
sponds to the particle being localized in one of the po- 
tential wells, to a flnitc voltage state, when the min- 
ima disappear and the particle runs down the poten- 
tial, at 7 ^ 1. If, like in our case, the temperature 
is finite, the particle may escape from the well already 
7 5- 1 , beein g th ern ially excited over the barrier 

H, MM, M MM, M m- 



The rate at which this process occurs depends on the 
barrier height 



ua = 2Ej 



■\/l^^7^ — 7 arccos(7) 



(4) 



and the small amplitude oscillation frequency of the par- 
ticle at the bottom of the well, 



Wo = 



= cjp (1 - 7 



(5) 



Here, cUp ~ y/2TTlc/^oC is the plasma frequency. For 
1 — 7 ^ 1 , Eq. ([H can be approximated as 



~ 77 4\/2 3/2 
wo ~ Ej—- (1 -7) ' 



(6) 



In the thermal regime the escape of the particle from the 
well occurs at a bias-current dependent rate of [2^, [13] 



T.. N ^0(7) 

r(7) = ^xp 



"0(7) 



(7) 



where a < 1 is a damping dependent prefactor. Knowing 
all relevant junction parameters and the bath tempera- 
ture T, one can calculate the escape rate F of a current- 
biased point-like J J. 



III. LONG JOSEPHSON JUNCTION WITH A 
FRACTIONAL VORTEX 

Here, the considered system, a long JJ with a frac- 
tional vortex, is more complex than a point-like J J. In 
order to account for the spatial phase variation one 
has to extend the picture of the particle in a tilted wash- 
board potential. First, to take into account the finite 
junction length, one replaces the point-like particle in a 
ID-potcntial by an clastic string in an extended poten- 
tial. Second, the k vortex in the LJJ corresponds to a 
kink in the string, as illustrated in Fig. [21 




FIG. 2: LJJ in the picture of a particle chain in an extended 
potential. 
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The depinning process of a fractional vortex at 7 < 
7c(/«) [see Eq. ([1])], described above, corresponds to the 
escape of this string out of a metastable well. To find 
the effective barrier height for this system, we map the 
multidimensional system to a point-like particle in an 
effective ID potential. By comparing the escape rate F, 
which we obtain by high resolution measurement of the 
depinning current, to the known escape rate of a point- 
like J J, one can make a statement about the effective 
barrier height. 

For an infinitely long J J, the effective potential for the 
escape of a fractional vortex close to the depinning cur- 
rent can be derived analytically, as shown in Ref. [3l|. 
Here, only the main results are summarized. The barrier 
height is given by 



with 



2 K 

- sin - (sin Pc - IcPc) 



3/2 



/3c = - - arcsin(7c). 



(8) 

(9) 
(10) 



and the depinning current 7c (k) given by Eq. ([T]) 
^ojcwXj /2t: = EjXj/L is the Joscphson coupling en- 
ergy per normalized unit length and jc the junction's 
critical current density assumed to be equal in and n 
parts. Here, we use and reserve uq as the barrier 
height for the homogeneous phase escape (short junction 
limit). The small oscillation frequency (O in this case is 
given by the eigenfrequency of the fractional vortex. In 
absence of a bias current the eigenfrequency is given by 

MM 




3cos2-). (11) 



For 7 7^ the analytical expression for luo(k,"/) is un- 
known, but can be approximated as 



^^o{i^,l) ~ cjo(k,0)(/ 1 
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(12) 



For LJJs of finite length one may use a numerical ap- 
proach and follow the ansatz described by Castellano et 
al. [13]. As in a point-like J J the thermal activation 
occurs over a barrier that is defined by the energy dif- 
ference between a minimum and the lowest adjacent sad- 
dle point/maximum of the potential surface. The corre- 
sponding phase distributions are solutions of the station- 
ary sine- Gordon equation 



fixx - sin [/i + kE{x)\ -1-7 = 0, 



(13) 



where Ji{x) is a step function that describes the position 
of the and k region 



for .T < 

1 for a; > 0. 



(14) 



For a linear junction geometry at zero applied magnetic 
field the phase obeys the boundary conditions 



(15) 



Here, /^(x) is the Josephson phase, which is a continuous 
function of cc, the spatial coordinate x is normalized to 
the Josephson penetration depth Aj and I = L/Xj is 
the normalized junction length. The subscript x denotes 
the derivate with respect to coordinate x. The potential 
energy of any state /i(a;) is given by 



U = E,\ 



1/2 
-111 



1 - cos[^ -I- kH{x)\ - 7Ai + ^ 



dx. 



(16) 



Let us denote the solutions corresponding to energy 
minima as fia{x) and the ones corresponding to sad- 
dle points/maxima as fisix). Then, the effective barrier 
height is 



u = U{pLs) ~ U{na), 



(17) 



which is defined for each pair of fig, fia- 

For an annular geometry it is more convenient to 
rewrite and solve the sine-Gordon equation p3p for the 
(discontinuous) phase 4>{x) = ii{x) + kH{x): 



sin (/} + "/ + kSx = 0, 



(18) 



where is the derivate of the (5-function. The vortex is 
placed at X = 0, which in annular geometry corresponds 
to a; = /, and the following boundary conditions are used 



m 

MO) 



0(0 + K, 
Ml)- 



(19) 



If K > 0, then the topological charge of the vortex is 
— K<0. Note, that all static solutions of Eqs. p8|) and 
also obey ^ 

7 = [(1 — cosk) cos 0(0) — sin K sin 0(0)] . (20) 

Here, 0(0) denotes the phase value at a; = 0. By inverting 
the signs of k and 0(0), Eq. (fT8|) and Eq. (|19p can also 
be applied to a vortex with positive charge. Eq. ([^0)1 is 
particularly useful, if a shooting algorithm is used to find 
solutions of Eqs. HH]) and P^ . 

For K = 0, the effect of junction length (and an ex- 
ternally applied magnetic field) on the phase escape in 
linear LJJs has already been discussed by other authors 
[27| . In the following we will summarize some of the re- 
sults, particularly important for the situation discussed 
in this paper. That is, we do not consider an applied mag- 
netic field. In short junctions, I < tt, the barrier height u 
is identical to uq for any 7, i.e., the value obtained from 
the short junction limit formula ([4]). The phase config- 
urations corresponding to minima and saddle points are 
homogeneous, i.e., fJ.{x) — const.. In longer junctions, 
I > TT, phase escape via additional, inhomogcneous saddle 
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solutions of lower potential barrier (jl7p results in w < uq. 
This is easily understood, as the energy required for a ho- 
mogeneous phase escape increases linearly with junction 
length, whereas it is constant [39| for the alternative pro- 
cess, that is the creation of a fluxon/antifluxon and its 
successive propagation along the junction. However, the 
nature and ultimate existence of the saddle solutions de- 
pends on bias current 7. As it turns out, at dc-bias close 
to the critical current, 7c — 7 ^ 1, which is of primary 
relevance in experiment, usually only one minimum and 
one saddle point solution remains. For junction lengths 
I < 9 and bias currents 7 ^ 0.997c, the barrier heights u 
even coincide with the short junction value uq. In annu- 
lar junctions the length at which the phase escapes via 
inhomogeneous phase configurations is 2tt instead of tt, 
as boundary conditions require the creation of a fluxon- 
antifl-uxon pair. As in linear junctions, the barrier height 
u coincides with uq for 7 ^ 0.997c and I < 12, which is 
the longest junction considered here. 

Now, how does the presence of a fractional vortex, i.e., 
\k\ > 0, affect the phase escape? In the following we first 
consider an annular geometry, where junction boundaries 
do not play a role. The influence of open boundaries on 
phase escape is discussed in SecEl Fig. [3] compares the 
effective barrier height u{k) at 7 = 0.997c(k) for different 
junction lengths. Here, u is normalized to E'j to avoid a 
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FIG. 3: Effective barrier height at 7 = 0.997c as a function 
of K for different junction lengtlis. The gray line corresponds 
to the analytic result from Eq. (|8]). 

normalization that depends on JJ length. The analytic 
result for Z = 00 is also depicted in the graph. As 

can be seen, the presence of a vortex facilitates the phase 
escape by lowering the effective barrier height. The larger 
the vortex, the easier is the escape. Fig. [3] also shows, 
that with increasing junction length, u[k) asymptotically 
approaches - the analytic result for an infinitely 

long junction. 

For further discussion it is important to note, that 
although u strongly depends on bias current, numeri- 
cal calculations of ^(7) have shown, that it scales as 
(1 — 7/7c)^''^ close to critical current with a pre-factor 
almost independent of 7. This is very convenient, as u, 



Uy [compare Eq. ^] and uq [compare Eq. ([6|)] exhibit 
the same asymptotic behavior for 7 — *■ 7c- With their 
relative proportions independent of 7, they can be easily 
compared with experimental data. 

IV. EXPERIMENTS 

For the experiments we used tunnel Nb-Al-AlO^^-Nb 
annular LJJs (ALJJ) with different junction radii R. The 
junction properties are specified in Tab. HI the specific 
capacitance of all junctions is C ~ 4.1/iF/cm^. The 
Josephson penetration depth Aj was estimated taking 
into account the idle region [3^. The normalized length 
of the ALJJ is given by I = 2ttR/\j. 
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72.6 
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8.6 
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TABLE L Junction parameters (annular) at T = 4.2 K. 




FIG. 4: Optical image (top view) of one of the investigated 
samples: ALJJ with two pairs of injectors. The right pair of 
current injectors was not used during experiment. 

To create an arbitrary k discontinuity in the junctions, 
a pair of tiny current injectors is used, as shown in Fig.lH 
The short section of the top electrode between points A 
and B of a length ^ Aj has an inductance ijnj. A cur- 
rent /inj passing through this inductance creates a phase 
drop K = iinj/inj27r/$o across the distance AB, i.e., the 
K discontinuity. To calibrate the injectors we have mea- 
sured the critical current Ic as a function of /inj for the 
left injector pair in Fig. H) The Ic (/inj) pattern is pre- 
sented for the shortest junction with i? = 30 /im in Fig. [5] 
and looks like a perfect Fraunhofer pattern in accordance 
to the theory [li, [13, [H] . The first minimum at I^^™ cor- 
responds to K = 2tt. The values /j^"" are given in Tab. [H 
Thus, for any /jnj, the corresponding value of k can be 
calculated as k = 27r/inj/ 

The escape of a fractional vortex in an ALJJ is ex- 
perimentally investigated by measuring the statistics of 
switching of the junction from the zero-voltage state to 
a finite voltage state. To do so, a bias current applied 



5 




0.02 



-10 10 

injector current / . (mA) 



FIG. 5: The dependence 7c(/inj) measured at T ~ 4.2 K (gray 
symbols) and the corresponding theoretical curve 7c (k) (con- 
tinuous line). 



to the junction is ramped up at a constant rate / and 
the current Ic at which the junction switches from its 
zero- voltage state is recorded [IHl- The probability dis- 
tribution P{I) of such switching currents is found by ac- 
cumulating a large number of measurements of Ic and 
generating a histogram. Using the obtained P{I) distri- 
bution, the bias current dependent escape rate can be 
reconstructed as [201 



r(/) 



AI 



ir+AiPindi'' 



(21) 



where AI is the width of the bins on the histogram. The 
shape of the histogram depends on the bias current ramp 
rate, the bath temperature T and damping [s^ . Higher 
T leads to a broader histogram shifted towards lower cur- 
rents. Note, electronic noise in the measurement setup 
has, to the first order, the same effect as the increase of 
T. Since the accuracy of our measurement is defined by 
the width of the escape histogram, special measures were 
taken to suppress electronic noise down to the level where 
the broadening caused by electronic noise is well below 
the one caused by the bath temperature. We mounted 
the sample in a copper box to shield it from electromag- 
netic radiation and in a cryopcrm cylinder to shield it 
from static magnetic fields. The bias lines contained 
three-stage low pass filters: two cold filters in close vicin- 
ity of the sample and one at room temperature. The 
measurement technique and setup are similar to those 
described in Ref. [11]. 

Following the analysis method used in Refs. [H, [13, 
[sH], we fit the escape rate r(/) calculated from the mea- 
sured histogram using Eq. (|2ip with the known escape 
rate of a point-like J J, see Eq. having the relative 
barrier height u/uq and 7c /cO as fitting parameters. The 
damping dependent prefactor was assumed to be a = 1 . 

The effective barrier height found in this way is shown 
in Fig. [6] for 3 different junction lengths at a bath tem- 
perature T = 4.2 K. Solid circles indicate experimental 
results and the black solid line corresponds to numerical 
simulations. For all junction lengths, the experimental 



0.00 



is 

3 0.01 



0.00 
0.02 



J! 0.01 - 



0.00 



- ' 1 ' 

/=5.0 

H — ^^ 






. 1=6.7 

H — ^^ 




1 ' - 


/=8.6 .-^ 

/ • < 

— t y % 


• 


1 ' 

••• 



-271 



-Iti Oti Iti 271 

discontinuity k 



FIG. 6; Effective barrier height u(k) for different junction 
lengths / at T = 4.2 K. Solid circles correspond to experimen- 
tal results, white lines to theoretical predictions for infinite 
0-K LJJ and solid lines to numeric simulations using the sine- 
Gordon equation in ALJJ of given length. Dashed lines show 
barrier heights for point-like junctions with rescaled 7c. 



data are in good agreement with simulations in particu- 
lar for large k values. In the range of small values of the 
discontinuity k, the deviation from simulations increases 
with junction length. There are two possible reasons ex- 
plaining the deviations. 

First, this might be due to junction inhomogencities. 
For large only the Aj-vicinity of the K-discontinuity is 
important for the escape, whereas for |k| — > the whole 
junction matters. The longer the junction is, the higher 
is the probability for inhomogencities, thus the deviations 
are more dominant for longer JJs. 

Second, our calculations are based on a single mode ap- 
proximation, which becomes less accurate and even fails 
in the limit |k| and / — > oo[3l[- This happens be- 
cause the eigenfrequency of the lowest eigenmode and 
the next ones tJi, 0^2 etc. are not well separated at the 
values of bias current where the escape takes place. At 
first, it seems that at 7 ^ 7c the eigenfrequency 0, 
while Li;i(7c) stays finite. Therefore there is always a re- 
gion around 7c where u;o ^ wi, regardless how small is 
the separation between and wi at low 7, say, at 7 = 0. 
However, the peak of the escape histogram (even in the 
case of quantum escape with the rightmost position of 
the histogram) is not sufficiently close to 7c in order to 
keep the condition loq <^ uji satisfied. Then it is impor- 
tant how nmch loq and lui arc separated at those values 
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of 7 that are most relevant for the escape process. This 
separation is roughly proportional to, but much smaller 
than 011(0) — a;o(0). In the limit / 00 and |k| 0, 
according to Eq. (jlip . wo(0) — *■ ujp, while the the next 
mode (the lowest mode in the plasma band) uji = LOp. 
Thus, the single mode approximation fails. To make it 
working again one has cither to increase |k| or to reduce 
I. By increasing |k| one shifts wo(k, 0) downwards from 
LUi = LOp, according to Eq. pT|) . By decreasing I one in- 
creases the wave vector fci = 27r/Z of the lowest mode in 
the plasma band, thus, shifting o>i(0) up as 



01(0) = UJp^l + kl = L0p^l + {2TT/IY. 

Thus, a single mode approximation starts working again 
for short (point-like) JJ and any k. Note, that participa- 
tion of higher modes in the escape makes the escape rate 
higher, i.e., the barrier height calculated using a single- 
mode formula lower, exactly as observed in Fig. [6l 

The corresponding barrier height of a point-like JJ is 
also depicted in Fig. [S] the dashed line indicates uq for 
point-like junctions with a critical current 7c(k)/cOj given 
by Eq. ([T|). If the only effect caused by the k vortex is 
the suppression of the critical current of the junction, the 
experimental data should coincide with the dashed line. 
However, the experimental data and simulations deviate 
from this curve, which becomes more pronounced with 
increasing junction length. This indicates a difference 
in the escape process of a point-like junction and a long 
junction with a k vortex. 

The white line in Fig. [S] corresponds to the analytic 
result from Sec. IIIII for infinitely long JJs. Comparing 
this result to simulations and experimental data, three 
regimes can be distinguished: a homogeneous regime 
(small light gray region in Fig. [5]), an intermediate 
regime (moderate gray region in Fig.lH]) and a vortex 
escape regime (large |«;|, dark gray region in Fig. [6]). In 
the homogeneous regime, the escape process corresponds 
to the one of a point-like junction: the phase string es- 
capes as a whole from the meta-stable state. The effective 
barrier height corresponds to uq in this regime. In the 
vortex escape regime, the escape-process is dominated by 
the vortex itself, and experimental data, simulations and 
theory coincide. 



height, but also with the junctions's resistance R assum- 
ing moderate to low damping [sj with 



(22) 



where Q = ojqRC is the quality factor. Using this 
method, the uncertainty of the experimentally obtained 
M(K)-dependence increases. Fig. [7] compares the barrier 
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FIG. 7: Effective barrier height m(k) for two different bath 
temperatures T in the thermal regime. Symbols indicate ex- 
perimental results, gray lines show theoretical predictions for 
the vortex escape and solid lines show numerical simulations. 

height determined at T = 500 mK for two annular JJs 
of lengths I = 6.7 and / = 8.6 with the already shown 
results at T = 4.2 K, c/. Fig. [B] The resistance values 
obtained from best fit are R = 480 and R = 433 51, 
respectively, corresponding to Q- values of ~2000 ... 3000. 
Note, the escape rate F only weakly depends on R, which 
leads to relatively large uncertainties in the determined R 
values (~ 50%). Though the method of determining u is 
less accurate for lower temperatures, the measured u{k)- 
dependence shows very good agreement with numerical 
simulations. 



Temperature dependence 

In the thermal regime, the effective barrier height u 
does not depend on the bath temperature T, see Eq. ([B]). 
For temperatures well above the crossover temperature 
to the quantum regime (23 |. which for our samples is 
~ 100 mK, the u(K)-dependencies measured at different 
temperatures should coincide. For lower bath temper- 
atures, T <C 4.2 K, the damping dependent prefactor a 
in Eq. ([7]) becomes important [35[. To account for this 
unknown factor, we fit the escape rate r(/) calculated 
from the measured histogram not only with the barrier 



V. LINEAR JUNCTIONS 

To study the effect caused by the junction geometry 
we have also measured linear LJJs with different lengths. 
As shown in Ref. [sst . the 7c (k) pattern of linear junc- 
tions differs from Eq. ([T]) due to the open ends. It is 
27r-periodic, has minima at k = tt -I- 27rn and depends 
on junction length. Calibration of the injectors follows 
the same method as described for annular junctions, how- 
ever the first minimum of the /c(/inj) pattern corresponds 
here to k = tt. Knowing the values Z™]", the correspond- 
ing K can be calculated as k = T^h-ai/ |-^iJJJ"| for any /inj. 
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Note, in contrast to annular junctions, 7c(k) depends on 
lengths I and has to be determined numerically. In addi- 
tion, only vortices with |k| < tt can be investigated with 
the experimental technique described above. For further 
details see Ref. j38| . 

For the experiments, we used Nb-Al-AlOx-Nb linear 
LJJs with different lengths, as shown in Fig. [51 The 
junction properties are listed in Tab. [TTl 



L 


w 




Aj 


/ 


rmin 
-'inj 


[pm] 


[fim] 


[A/cm^] 


[fim] 




[mA] 


120 


5 


84.2 


53.5 


2.3 


4.9 


240 


5 


78.5 


55.4 


4.4 


5.0 



TABLE II: Junction parameters (linear) at T = 4.2 K. 




FIG. 8: Optical image of one of the investigated LJJs with 
one pair of injectors. 




I I I I I I I I I 

-LOti -0.5% O.On O.Sn 1.0% 

discontinuity k 



FIG. 9: Effective barrier height as a function of k for linear 
junctions with different lengths at T = 4.2 K. Solid circles 
show experimental results, gray lines show the theoretical 
prediction and solid lines show numerical simulations. The 
dashed line indicates the barrier height of an annular junc- 
tion of the same junction length. 

In Fig. [9] the results for linear LJJs are shown for 
T = 4.2 K. Experimental data and numerical simulations 
show a very good agreement. On the other hand, a com- 
parison of experimental data and theory shows, that even 



in the case of the longest junction with I = 4.4 the ex- 
perimental barrier height is considerably lower than the 
asymptotic value u„(k;) even for high values of k. Note 
that for annular junctions all three (experimental data, 
numerical simulations and analytic calculations) are in 
good agreement. This discrepancy is due to the junc- 
tion geometry: the escape in linear junctions is usually 
dominated by the nucleation of a fluxon-antifluxon pair 
at the edges. Only at very large / a crossover to vortex 
activation takes place. At the same time, the observed 
K-dcpcndcncc suggests that boundary dynamics are af- 
fected by the presence of a fractional vortex. This be- 
comes clear, if one compares numerically simulated u{k)- 
dependencies for two junctions with different geometry, 
as shown for the longest junction with I = 4.4: the dashed 
line in Fig. [5] corresponds to an annular junction of the 
same length. The obtained barrier heights differ in shape 
and absolute value. Fig. [10] compares simulated barrier 



1 ' 


1 « 




7 A ■ 




/ J'K/ 































-LOti -O.Sti O.Oti 0.5ji LOti 

discontinuity k 

FIG. 10: Numerically calculated effective barrier heights at 
7 = 0.997c as a function of k for linear junctions with different 
lengths. The gray line corresponds to analytical results. 

height dependencies on k of linear junctions with different 
lengths to the theoretical prediction u^(k) for infinitely 
long JJs. Only in rather long junctions [l > 9) and for 
high values of k the escape mechanism caused by the frac- 
tional vortex becomes important. For I < 9 and small k 
the edge effects doininate the phase escape. 



VI. CONCLUSION 

We have experimentally investigated the thermal es- 
cape of an arbitrary fractional Joscphson vortex close to 
its depinning current 7c (k) in long JJs of different geom- 
etry. A comparison of our results with numerical calcula- 
tions based on the perturbed sine-Gordon equation shows 
very good agreement. Depending on the junction's geom- 
etry, our results coincide with theoretical predictions in 
different regimes of k, except in the limit k — > 0, Z — > oo, 
where our model, which is based on a single mode ap- 
proximation, fails [3]| . 
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